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Abstract 

This  report  includes  the  calibration  results  of  the  Lab(k)-PTB  TWSTFT  link  and  closure  measurements  of  the  BIPM-Lab(k)- 
BIPM  tour.  During  10-20  Feb.,  2015  (DOY  41-51,  MJD  57063-57073),  the  BIPM  Standard  Travelling  Calibration  Station 
(Stde)  visited  Lab(k)  in  order  to  calibrate  the  Lab(k)-PTB  TWSTFT  link  for  UTC  generation.  This  work  follows  the  TWSTFT 
Calibration  Guideline  for  UTC  Time  Links  [1], 


The  staff  and  the  calibration  setup  at  Lab(£) 


Annex  II  is  neither  mandated  by  the  guideline  [1]  nor  a  part  of  it  but  gives  an  example  how  to  perform  a 
TWSTFT  link  calibration  with  a  GPS  calibrator  and  prepare  the  report.  Acronyms  (BIPM,  USNO  etc.), 
coordinator,  equipment,  images,  figures,  measurement  data  etc.  in  this  document  are  all  fictional. 
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Notation 

UTCp:  the  UTC(£)  point  at  Lab(£).  Here  after  the  k  stands  for  Lab(k),  the  laboratory  to  be  calibrated 

Link:  a  time  link  is  a  clock  comparison  result  using  a  particular  technique,  e.g.,  a  link  of  GPS  C/A,  P3,  PPP  or  GLONASS  or 
TWSTFT  or  TWOTFT.  A  UTC  link  at  present  is  a  time  link  between  Lab(£)  and  PTB 

StdB:  The  GPS  travelling  calibrator.  In  this  Lab(k)  calibration  tour,  the  BIPM  standard  traveling  calibration  station 

(calibrator)  consisting  of  N  (>2)  GNSS  receivers+antennas+cables  and  PPS/frequency-distributors.  It  is  a  pre-cabled 
black  box  calibrator  with  unknown  but  constant  total  delay  during  a  calibration  tour 

Total  Delay:  The  total  electrical  delay  from  the  antenna  phase  center  to  the  UTCp  including  all  the  devices/cables  that  the 
satellite  and  clock  signals  pass  through.  It  numerically  equals  the  sum  of  all  the  sub-delays.  The  total  delay 
uncertainty  determines  the  UTC  time  transfer  uncertainty 

METODE2:  MEasurement  of  TOtal  DElay,  the  BIPM  calibration  system  composed  of  related  methods  and  equipment  (StdB) 
for  the  generation  of  UTC-UTC(L)  in  Circular  T  [1] 

CM:  The  METODE  total  delay  correction.  It  should  be  subtracted  from  the  GPS  data,  e.g.  RefGPS-CM  in  CGGTTS,  -CM  in 
Clb_GNSS.Lst  file;  and  added  to  the  CALR  of  the  ITU  TWSTFT  data  of  the  Lab(C)  side.  Because  the  PTB  is  taken 
as  the  reference  of  the  calibration,  a  GNSS  time  link  correction  is  equal  to  the  classic  GNSS  equipment  calibration 
correction  [8] 

uA,  uB:  type  A  and  type  B  uncertainties  (l-o) 

Um:  Total  uncertainty  of  the  total  delay  correction  Cm; 

CCD:  a  difference  of  two  system’s  data  that  have  a  common  clock 

DCD:  double  difference  of  two  independent  measurements  of  clock  differences,  for  example  via  TWSTFT  and  GPS. 

Tour:  a  calibration  tour  is  a  round  trip  calibration  campaign  with  start  and  closure  measurements.  It  may  include  several 
laboratories 


2  METODE  was  proposed  in  the  frame  of  the  BIPM  pilot  project  (2011-2014)  aiming  at  unifying  the  UTC  time  link  calibration  with  an 
uncertainty  <  2  ns  [2-4],  It  is  composed  of  a  time  link  calibration  scheme  with  the  calibrator  denoted  StdB. 

This  document  describes  a  typical  TWSTFT  link  calibration.  If  we  replace  the  TWSTFT  link  by  a  GNSS  link  or  a  optical  fiber  (OF),  it 
becomes  a  GNSS  or  an  OF  time  link  calibration.  This  calibration  becomes  a  classic  receiver  calibration  if  the  link  includes  the  UTC  network 
pivot  (PTB),  whose  absolute  calibration  error  is  assumed  to  be  zero.  The  calibration  correction  Cm  can  be  converted  to  classical  corrections 
of  the  Internal  Delay,  INTDLY(L1/L2),  by  removing  the  CABDLY  and  REFDLY.  This  introduces  extra  uncertainties.  In  consequence,  the 
uncertainties  of  the  INTDLY(L1/L2)  maybe  larger  than  3  ns  [3,15], 
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1  Summary 

According  to  the  TWSTFT  guideline  [1],  a  TWSTFT  link  calibration  campaign  is  carried  out  using  a  mobile 
TWSTFT  ground  station  or/and  a  GPS  calibrator  that  are  circulated  among  several  time  laboratories  contributing 
to  UTC.  This  report  confines  itself  to  the  specific  measurement  of  Lab(k)-PTB.  A  similar  calibration  tour  of 
NIST-PTB  has  been  made,  the  consistency  of  that  report  with  previous  TWSTFT  calibrations  of  the  Lab(k)-PTB 
and  Lab(k)-NIST  links  will  be  presented  elsewhere,  as  will  be  the  issues  involved  with  the  application  of  these 
calibrations  to  the  GPS  receivers  at  the  two  sites. 

1.1  General 

This  report  includes  the  calibration  results  of  the  Lab(k)-PTB  TWSTFT  link  and  closure  measurement  of  the 
BIPM-Lab(k)-BIPM  tour  with  the  BIPM  standard  travelling  calibration  station  (StdB).  During  10-20  Feb.,  2015 
(DOY  41-51,  MJD  57063-57073),  the  StdB  was  installed  at  Lab(k),  United  States  Naval  Observatory,  3450 
Massachusetts  Ave,  NW,  Washington  DC  20392  USA.  The  goal  was  to  calibrate  the  Lab(k)-PTB  TWSTFT  link 
for  UTC  generation.  This  work  and  this  report  follow  the  TWSTFT  Calibration  Guideline  for  UTC  Time  Links 
[1]. 

As  part  of  the  BIPM  Pilot  Project,  the  METODE  was  developed  to  unify  the  UTC  time  link  calibrations  with  a 
calibration  uncertainty  uB  <  2  ns  [2-5],  Since  2013,  the  StdB  has  visited  the  UTC  labs  OP,  PTB,  PL,  AOS,  TL, 
NMIJ,  NICT,  NIM  (BSNC),  and  ROA  [6];  experiments  were  made  also  at  the  BIPM,  NIST  and  Lab(k)  [7,8]. 
The  two  StdB’s  visits  to  PTB  in  June  2013  and  Aug.  2014  [6]  allow  transferring  the  calibration  of  the  PTB 
master  receiver  to  the  Lab(£).  The  difference  of  the  two  visits  is  0.03  ns.  This  and  the  closure  measurements  at 
BIPM  prove  the  long-term  stability  of  the  StdB. 

The  requirements  for  the  setup  and  computations  can  be  found  in  the  BIPM  guideline  [5].  Taking  into  account  of 
the  starting  and  closure  measurements  at  the  BIPM,  we  compute  the  calibration  corrections  for  the  UTC 
TWSTFT  time  links  between  Lab(k)-PTB.  Since  this  link  has  been  recently  calibrated  with  TWSTFT,  this 
supplies  a  supplementary  and  official  GPS  time  link  calibration. 

GPSPPP  solutions  are  used  for  this  calibration,  and  hereafter  the  UCT(£)  or  Lab(/t),  refers  only  to  k=  Lab(k). 

1.2  Summary  of  the  main  result 

Table  1  displays  the  calibration  for  the  time  link  corrections  (CM)  for  the  TWSTFT  links  on  the  baseline  Lab(k)- 
PTB,  see  Table  1. 

Table  1  The  total  delay  correction  for  the  TWSTFT  time  link  Lab(k)-PTB 


Lab 

Time  Rcv/Link 

CmAis 

Um 

ITU  CI 

S 

Lab(k) 

TWSTFT:  Lab(k)-PTB 

+0  9* 

1.5  ns 

888-2015 

1 

*  In  the  files  TWLABK57.070,  and  TWPTB57.070,  we  have  the  corresponding  CALR=-488.0  ns  for  Lab(k) 
and  +488.0  ns  for  PTB.  The  ESDVARs=  are  kept  zero  and  unchanged,  cf.  Section  4.2. 

Figure  1  shows  the  data  of  the  related  links.  Of  the  available  GPSPPP  links,  only  USN6-PTBB  is  illustrated  here. 
As  shown  in  the  plot,  the  Lab(k)-PTB  TWSTFT  and  the  GPSPPP  links  are  close  to  each  other  but  0.9  ns  and  1.5 
ns  lower  than  the  StdB-PTBB  link. 


violet:  PPP  Lab(k)-PTB;  Triangle:  TWSTFT  Lab(k)-PTB;  Red:  PPP  StdB-PTBB  (StdB  is  mean  of  BPOU/blue  and  BPIC/black) 
Figure  1  The  time  links  on  the  UTC  baseline  Lab(k)-PTB  during  the  calibration  period 
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1.3  Th  ?  Combined  Uncerta  nty 

The  total  uncertainty  (Um)  of  the  C  [  is  composed  of  [3-4]: 

•  PPP  Measurement  uncertainty  (ua)  of  StdB-UTC(k):  0.1  "0.3  ns: 

•  PPP  Measurement  uncertainty  (ua)  of  UTC(k)-  JTC(PTB):  0. 1-0.3  ns 

•  TWST  FT  Measurement  uncertainty  (ua)  of  UTC(k)-UTCfPTB):  0.2-0. 5 

•  Instabi  ity  and  the  sub-delay  measurement  uncertainty  of  the  reference  at  .ab(k):  0.5-0. 7  ns 

•  Instabi  ity  of  die  traveling  receivers:  0.5-1 .0  ns; 

•  Others:  0.3-0. 6  ns 

Hie  UM  as  estimated  from  the  root  sum  square  (RSS)  of  the  e  errors  is  hence  (0.8-1. 5)  ns  (la).  However,  the 
unknow  .  and  systematic  errors  are  likely  to  be  higher  dian  us  lal  and  therefore  we  take  our  uncertainty  to  be  1 .5 


If  only  one  GPS  receiver  calibrato  •  component  of  the  MET  'DE  is  used,  the  instability  would  be  factor  of  v2 
higher,  a  id  //B  =  s/2  xl.5  ns  =  2.1  ns. 

Other  independent  studies  [6-15]  prove  that  th  •  calibration  mcertainty  of  1.5  ns  or  even  below  1  ns  [16]  is 
attainabl ;.  Here  each  system  can  perform  a  calibration  without  sharing  any  common  part  with  the  other.  It  is  best 
to  have  at  least  two  receivers  of  different  types.  Hiis  may  inc  ease  the  measurement  discrepancies  bu:  improves 
the  uncertainly  computation  as  well  as  die  robustness  of  the  calibration  result. 


2  Sta  idard  setups  of  t  le  StdB  during  a  calibrati  n  tour 


The  setu )  of  die  StdB  is  shown  in  the  Figure  2  [5].  The  cable  €166  was  dire  :tly  connected  to  the  UTC(b). 


Lab(/fc) 

BIPM 

BIPM 

ante  ina(s) 

antenna  1 

antenna  2 

Lab(A')  clock 


10  MHz  distributor 


PPS  distributor 


UTC(*)L  - 


Time  transfer 
receiver  of  Lab(£) 


Lab(k) 

equipment 


lOMhz 

C149 


1PPS 

066 


10  Ml  Iz  distributor 


ITS  ili'lrihtuor 
BPIC/?olarX3 


BPOU/GTR50 


TIC 


Fi  jure  2  Setup  if  the  BIPM  StdB  at  a  UTC  ( k  ) 

(Tlie  BIPM  devices  incl  ding  cables  are  shown  in  blue.  Lab(A)’s  equipment  are  shown  in  black) 

By  the  d  :finition  of  the  METODE  JTC  time  lin  :  calibration  correction  [2,3].  we  have  the  following  steps: 

•  We  s  art  fr  om  BIPM; 

•  We  et  the  PTB's  master  GP>  receiver  (PTBB)  as  the  reference  of  the  calibration  and  its  :alibration 

cor  ection  to  be  zero: 

•  We  align  the  StdB  to  PTBB.  i.e.  the  BP0U  an  1  BP1C  in  St  B  are  to  be  corrected  -5.2  ns  and  -3.6  ns  [3]: 

•  The  >tdB  goes  to  the  Lab(A').  and  makes  m  tasurements  side  by  side  widi  the  TW  >TFT  ground  station  of 

Lab(A').  They  use  the  same  reference  signals  of  UTC(A); 

•  We  make  the  closure  measurem  int  at  BIPM: 
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•  We  compute  the  double  clock  differences  (DCDs)  as  shown  in  equation  (1).  Each  data  point  is  the  result  of 
the  difference  of  a  TWSTFT  value  and  the  interpolation  of  the  2  adjacent  PPP  values  (computed  every  5 
min)  or  P3  values  (computed  about  every  16  min).  P3  technique  is  not  the  best  option  to  cany  out  DCD, 
even  worse  for  long  baseline,  nevertheless  the  differences  with  respect  to  PPP  results  are  normally  below 
0.5  ns.  The  corresponding  equation  for  the  DCD  is: 

Cm=DC'D  =Liuk(PPP)-Link(TW)  =[UTCppp(A-)s<dB-UTCppp(PTB)]  -[UTC(£)-UTC(PTB)]tw  (1) 

here  the  GPSPPP  data  in  the  first  bracket  are  taken  while  StdB  is  at  site  k.  The  UTCxw(^)  is  measured  by  die 
TWSTFT  equipment;  The  no-zero  DCD  is  the  calibration  correction  to  the  link  Lab(A')-PTB. 

To  average  out  the  diurnal  effects  and  measurement  noise,  5-7  days  of  continuous  measurements  is  requir  ed. 


3  Setups  at  the  Lab [k] 

The  setup  and  the  1-PPS  IN/OUT  measurements  at  Lab(k)  are  illustrated  in  the  Figures  3.1  and  3.2.  See  also  the 
photos  on  the  cover  page.  The  RefDly  determination  is  critical,  and  is  the  only  value  that  must  be  meastued  in 
both  laboratories.  Although  not  difficult  in  principle,  subtle  impedance  matching  issues,  reflections,  and  even  the 
choicee  of  measurement  technique  could  can  affect  the  measurement  [2,15].  In  the  BIPM  StdB.  a  time  interval 
counter  (TIC)  is  shipped  with  to  reduce  the  impact  of  the  bias  in  the  sub-delay  measurements  (Figure  2).  In  the 
setup  of  Figure  3.1.  the  RefDly  of  the  StdB  is  66.2  ns. 


BIPM  calibration  unit 


Figure  3.1  Setup  of  the  StdB  at  the  Lab(k)  T/F  laboratory 


|  66.2ns  | 


USNO  MC 
1PPS 

1  PPS  Distributor 

• - Cable 

1  (Cl  57) - % 

BNC  adaptor 

Cable  3<C1«3) 


Null  measurement  f51 

| - Cable 

STARti  1  STOP 

Time  Interval  Counter 

Figure  3.2  BIPM  BP1C  PPS  IN/OUT  measurements 
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Figure  3.2  shows  the  BIPM  BP1C  PPS  IN/OUT  measurement  on  10/2/2015  before  the  calibration  measurement 
started.  On  20/2/2015.  another  measurement  wras  made  after  the  measurement.  The  difference  was  0.05  ns  and  is 
negligible. 


Table  3.1  is  a  summary  of  the  receiver  and  the  antenna  information  directly  used  in  the  calibration  data 
processing.  Table  3.2  lists  the  present  sub-delays  before  the  calibration.  They  will  be  used  as  the  starting  values 
for  the  calibration  computations. 


Table  3.1  Hie  receiver  and  antenna  information 


No. 

Receivei 

Type 

Antenna 

Antenna  code 

Note 

1 

BP0U 

GTR50 

NOV702GG 

NAE07 190046 

BIPM  StdB 

2 

BP1C 

Sept.  Polarx3 

ASH701945E  M 

2000785 

BIPM  StdB 

3 

PTBB 

Ashtech  Z12T 

ASH700936E  SNOW 

CR15930 

Master 

# 


Table  3.2  The  sub-delay  information  (in  CGGTTS  header  etc.)/ns 


No. 

Receiver 

IntDly(Ll) 

IntDly(L2) 

D 

f 

CabDly 

RefDly 

Co* 

Cl* 

C2* 

rrn 

Mil 

TotalDly 

Note 

1 

BP0U 

-66.2 

-20.8 

5.2 

-81.8 

2 

BP1C 

-66.2 

225.2 

3.6 

-203.3 

-6.5 

-412 

before 

3 

PTBB 

304.5 

318.9 

282.252 

301.7 

75.3 

508.65 

*  Co,  Cl,  C2,  C3  are  the  sub-delays/coirections.  We  use  only  the  Total  Delay  for  the  link  calibration 


4  Data  reduction  and  analysis 

We  use  the  equation  (1)  to  compute  the  total  delay  calibration  correction  Cm  through  the  DCD  of  the  TWSTFT 
and  the  GPSPPP  links. 

4.1  GPSPPP  solution 

The  RINEX  files  (including  PTBB)  were  edited  and  corrected  for  cycle  slipsith  the  program  Teqc  before  the  PPP 
processing.  For  the  Novatel  receiver,  the  bias  Cl-Pl  was  taken  into  account  using  CC2nonCC. 

Figures  4.1.1a  and  4.1.1b  show  the  DCD  of  the  BP0U  and  BP1C  and  the  TDev.  The  measurements  were 
disturbed  by  unknown  reasons.  The  DCD  scatter  up  to  1  ns.  Although  most  of  the  deviations  should  be  averaged 
out,  the  possibility  of  a  systematic  bias  in  one  of  the  two  receivers  cannot  be  ruled  out  by  divided  by  2.  The  mean 
value  is  0.357±0.200  ns. 

The  red  triangles  in  the  following  figures  are  the  day-averaged  values. 


_ , _ , _ , _ +  , _ , _ , _ ’  ,  _ , _ 

57073  57074  57075  57076  57077  57078  57079  57080 
MJD 

Figure  4.1.1a  DCD  of  BP1C-BP0U.  Av  =  0.357±0.200  ns 
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2.0  2.5  3.0  3.5  4.0  4.5  5.0  5.5 

Averaging  time 

Figure  4.1.1b  Tdev  of  the  DCD  in  Figure  4.1.1a 


4.3  The  calibration  of  the  TWSTFT  link 

The  raw  data  of  the  GPS  and  TWSTFT  between  MJD  57075-57081  were  used. 

Figures  4.1.2a  and  4.1.2b  depict  the  CCD  and  DCD  of  the  TWSTFT  minus  GPSPPP  links  over  the  baseline 
Lab(k)-PTB.  Here  and  below,  the  black  cross  is  TW  link  and  blue  circle  the  PPP  link.  The  DCD.  i.e..  the 
calibration  corrections,  are  0.79  ns  and  1.07  ns  as  measured  by  the  BP0U  and  BP1C  respectively.  Their  average 
TWSTFT  link  calibration  correction  is  Cm=  0.93  ns  wTith  an  uncertainty  of  1.5  ns.  If  applied,  it  should  be  added 
to  the  CALR  in  the  ITU  file  on  the  Lab(k)  side  and  subtracted  from  the  PTB  ITU  file.. 


25  T 


57075  57076  57077  57078  57079  57080 

MJD 
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Figure  4.1.2a  DCD  of  TWSTFT  and  PPP  (BPOU-PTBB)  links  of  Lab(k)-PTB.  Av=0.79±0.309  ns 


57076.0  57077.0  57078.0 

MJD 


57079.0 


57080.0 


Figure  4.1.2b  DCD  of  TWSTFT  and  PPP  (5P7C-PTBB)  links  of  Lab(k)-PTB.  Av=1.07±0.339  ns 


Tlie  present  CALR=  -488.884  ns  with  ESDVAR=0.000±0.0  ns  in  the  ITU  file  TWLABK57.070.  We  have  then 
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the  METODEcalibrated  CALR=  -488.88+(0.93)=-487.95  ns  with  present  ESDVAR=0  kept  imchanged  in  both 
sides  of  PTB  and  Lab(k). 

Tliis  coiTection  should  be  subtracted  from  the  ITU  TWSTFT  data  format  file  of  the  PTB  side  but  added  to  that 
of  Lab(k)  side.  The  Job  of  the  BIPM  Tsoft  Menu  Y20  for  this  calibration  correction  (active  Calib)  is: 

Calib. _ :  S=0  CALR=  00000.000  ESDVAR=  00000.000  !  CALR=-488 . 0=-488 . 9+0 . 9/ITU  CI=394 

Calib.  PTB 01  OSNOOl :  S=0  CALR=  +488.0  ESDVAR=  00000.000  'subtracted  from  ITU  TWPTBmj.ddd 

Calib.  USNO01  PTB01 :  S=0  CALR=  -488.0  ESDVAR=  00000.000  ladded  to  ITU  TWUSNOm j . ddd  files 

Note  here  that,  usually  the  ESDVAR  should  be  set  to  zero  after  calibration. 


4.2  The  TWSTFT  and  GPSPPP  links  after  the  calibrations 


The  calibration  of  the  GNSS  time  transfer  facility  is  not  the  goal  of  this  calibration  tour.  We  give  the  following 
GPSPPP  and  TWSTFT  time  link  comparison  result  only  as  complementary  information. 

Figures  4.2.  la/b/c  show'  the  TWSTFT  and  USN6  GPSPPP  links  after  that  the  new  CALR(Labk)=  -488.0  ns  (ITU 
0=394/  S=l)  and  the  new  INTDLY(L3)  =  +1.5  ns  are  applied,  cf.  Table  1  and  [8].  The  mean  of  the  differences 
is  0.044  ns  ±0.284  ns.  Diumals  in  both  GPSPPP  and  TWSTFT  present  as  shown  in  the  Figure  4.2.1a. 


2.8- 

2.6- 

2.4;. 

2  2"  ‘ 
c/>  . 

-5  2.0- 

8 1 8  + 

O  1.6 

1.4 

1.2- 

1.0- 


Fig.  4.2.1a  (72  Ep) 

TW  link  (black  cross)  and  GPSPPP  link  (blue  circle) 


57076.0 


57077.0 


57078.0 

MJD 


57079.0 


57080.0 


u 
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-0 

-0 


57076  57077  57078  57079  57080 

MJD 


Figure  4.2.1  DCD  of  Lab(k)-PTB  link  comparison.  Both  TW  and  PPP  links  are  calibrated 


Fig.  4.2.1b  TW  link  -  GPSPPP  link 
Mean=0.044  ns  ±0.284  ns 
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5  Stability  of  Stcb  and  closure  at  BIPM  before/after  the  Lab (k)  tour 

The  final  computation  should  be  made  after  the  closure  measurement  which  controls  the  stability  of  the  StdB. 
The  Table  5.1  gives  the  GPSPPP  closures  at  BIPM  before  and  after  the  visit  to  Lab(k).  The  stationary  BIPM 
receiver  GTR50  BPOT  is  taken  as  a  reference.  On  average,  the  closure  of  the  StdB  is  -0.2  ns  for  the  two  receivers 
and  is  negligible'.  The  two  travelling  receivers  in  the  StdB  are  separated  each  other  of  0.4  ns.  The  StdB  is  hens 
stable  during  the  calibration  tour. 

The  instability  of  the  StdB  is  no  bigger  than  0.5  ns  since  its  last  visit  to  PTB  in  Aug.  2014. 

Table  5.1  The  PPP  closures  at  BIPM  before  and  after  the  visits  to  Lab(k)  vs.  the  GTR50  BPOT 


Period 

BPOU-BPOT/ns 

BPIC-BPOT/ns 

BPIC-BPOU/ns 

Mean  vs.  BPOT/ns 

57050-57056 

-0.4±0.2 

0.0±0.2 

0.4±0.2 

-0.2±0.2 

57090-57096 

-0.2±0.2 

0.2±0.2 

0.4±0.2 

-0.0±0.2 

Old-New  closure 

-0.2 

-0.2 

0.0 

-0.2 
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